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Abstract 

We present a semi-analytic method to calculate the dispersion curves and the group velocity of photonic 
crystal waveguide modes in two-dimensional geometries. We model the waveguide as a homogenous strip, 
surrounded by photonic crystal acting as diffracting mirrors. Following conventional guided-wave optics, 
the properties of the photonic crystal waveguide may be calculated from the phase upon propagation over 
the strip and the phase upon reflection. The cases of interest require a theory including the specular order 
and one other diffracted reflected order. The computational advantages let us scan a large parameter space, 
allowing us to find novel types of solutions. Keywords: photonic crystal waveguide; slow light; tailored 
dispersion; dispersion engineering 



1. Introduction 

Slow light has a number of important applications PQ [21 [3] , many of which stem from the strong electric 
field which is associated with a particular energy flow. A key geometry for the generation and observation 
of slow light are photonic crystal waveguides (PCWs) in which the dispersion relation, and hence the group 
velocity, can be tuned very precisely. 

One of the challenges in this field is the design of a structure with desired dispersion characteristics, 
typically involving both the group velocity and the dispersion. These can be engineered by varying the 
parameters of the first rows of the photonic crystal (PC) closest to the waveguide, for example the holes' 
radii or their transverse or lateral positions @]. In this way, waveguide modes with low group velocity 
which weakly depend on frequency have been found [5] |6]. The actual design procedure is often based 
on a combination of experience and fully numerical optimization. This optimization typically starts with a 
two-dimensional geometry, which is followed by full three-dimensional slab geometry simulations in the most 
promising part of parameter space. The properties of the three-dimensional geometry, such as bandwidth 
and group index, are well predicted by those of the two-dimensional geometry [5 . 

Here we consider an efficient and conceptually convenient method for the calculation of the dispersion 
relation of waveguides in two-dimensional PCs. The key to our approach is that the PCW is considered 
to be a uniform strip which is surrounded by two mirrors, similar to a conventional planar waveguide as 
illustrated in Fig. [T] The key difference, of course, is that in planar waveguides the mirrors rely on total 
internal refiection, whereas in a PCW they rely on a combination of Bragg reflection from the PC mirrors 
and total internal reflection from the PCW's faces. Nonetheless, as we show here, the type of methods 
developed to analyze planar waveguides can be brought to bear on the calculation of the properties of 
PCWs. This leads to a numerical method which is efficient in that the width of the waveguide may be 
varied with little computational cost, and which provides some physical insight. The design method for 
planar waveguides has been applied to PCWs previously. Sauvan and Lalanne used it to explain the PCW 
loss due to a finite cladding [7] . Istrate and Sargent concentrated on PCW-based devices such as directional 
couplers and cavities [8]. Chen et. al. used it to determine the modes of a waveguide with a photonic crystal 
cladding and an air cover [9]. None of these authors carried out a type of systematic slow-light analysis 
reported here. 
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Figure 1: Schematics of the two-dimensional geometries discussed, (a) Conventional planar geometry in which the mirror 
reflectivity stems from total internal reflection, (b) Scalar treatment of PCW in which the PC mirrors only reflect the specular 
order, (c) Vector treatment of PCW, in which the PC mirrors reflect the specular order and one diffracted order. 



In Sec. [2] we formulate the theory. In Sec. |2.1| we start with the case where only the specularly reflected 
order propagates, which applies to low-refractive index PCWs. As the PC mirrors serve as gratings, we 
formulate in Sec. |2.2| the theory when the PC mirrors reflect two propagating orders, which can be applied 
to practical PCWs. Some of the results for this are derived in the Appendix. In Sec. [3] we demonstrate the 
utility of our method and compare it to full three-dimensional numerical simulations. In Sec. [4] we discuss 
and conclude our findings. 

2. Description of the method 

2.1. Low-refractive index waveguides-scalar case 

It is well known that the resonance condition for a mode in a planar waveguide (see Fig. [ija)) is given 

by m 

2nkd COS -d + 2ip = 2rmr , (1) 

where d is the width of the waveguide, n is its refractive index, and fc = w/c is the wave number of the light 
in vacuum, with uj the frequency and c the speed of light in vacuum. Further, is the angle of the wave 
with the normal of the mirror and </? is the phase shift upon reflection (i.e., r = exp(i(p)). Equation ([T]) 
expresses the fact that the existence of a mode requires that the field is unchanged when traveling back 
and forth between the mirrors, which includes the phase accumulated due to the propagation and that 
due to two reflections off the mirrors, corresponding to the first and second terms on the left-hand-side of 
Eq. ([I]), respectively (see also the Appendix). The standard treatment for planar structures then continues 
by substituting the expression for ip = (fid) for total reflection which follows from the Fresnel coefficients, to 
find a transcendental equation for the waveguide modes. In contrast, in our treatment for PCWs (Fig.[l];b)) 
we obtain ip = ip{'d,ij) from a numerical calculation. We are also not only interested in the dispersion 
relation, but also in the group velocity. To this end we note that the propagation constant /3 = nfcsin??. 
The group velocity can thus be written as 

Ua = — = n ( sin-d + cost? I fc-rr I I . (2) 
«9 V V dkJJ 

We now implicitly differentiate Eq. ([T]) with respect to k to find the expression for kd'd/dk 

^dd _ ndcost9+|^|^ _ dtcost? , . 



The key elements, and the only elements which depend explicitly on the properties of the mirrors, are the 
partial derivatives of the phase upon reflection with respect to frequency and incident angle. The phase 
shift for reflection off an ideal metal is tt, and thus dt = dw = nd for a metal waveguide. In the most general 
situation the parameters d^ and d^. can be interpreted as follows: the partial derivative dip/d-d enters the 
expression for the Goos-Hanchcn shift, a sideway shift if the incident beam has a finite width, and which 
can be thought of as a correction to the position of a mirror lOj. Thus, d^ is the waveguide width corrected 
for the Goos-Hanchen shift. The inclusion of the factor n indicates that d^, corresponds to an optical path 
length. Similarly, the partial derivative dip/dcu is associated with the group delay, an apparent shift to the 
position of the mirror stemming from finite pulse lengths and dr is the waveguide width corrected for 
this effect. 

From Eqs ^ and ([3| we see that Vg ^ when the denominator in ([3| vanishes. Since the first term 
in the denominator in Isl) is positive, this means that dLp/dd needs to be small and positive. In regular 
dielectric waveguides dip/d-d is negative and so the group velocity does not vanish. In PCWs, Vg = both 
at the edge and at the middle of the Brillouin zone (BZ). We have found the general conditions for which 
Vg — Ohy considering the properties of the mirror only, but we do not consider these here, since, as discussed 
in the next paragraph, they apply only to waveguides made of low-index material. 

Results ([2]) and ([s]) are interesting but we have implicitly assumed the wavenumber in the waveguide to 
be small. This arises from the assumption that an incoming plane wave leads to a single reflected plane 
wave as in Fig. [ijb). However, each of the sides of the waveguide acts as a reflection grating, which, at 
sufficiently high frequencies, leads to diffracted orders, the angles of which are given by the grating equation 
[S]. For PC bandgap mirrors, result ([3| is valid only for an air waveguide, or waveguide made of another 
low-index material, sandwiched between the two PC mirrors. For the key case in which the waveguide is 
made of silicon or another high-index material, the frequencies for which Eq. ^ holds correspond to a band, 
rather than to the PCs' photonic band gap. To treat the more interesting case it is sufhcient to consider 
frequencies at which the mirror has one diffracted order in addition to the specularly reflected order, as 
discussed in Sec. 12.21 



2.2. High-refractive index waveguides-vector case 

Here we consider the case in which the mirrors have two reflected orders (labeled and —1), the specular 
order and one other (see Fig. [ijc)). In this situation the reflection coefficient is replaced by a 2 x 2 reflection 
matrix, in which the diagonal elements give the specular reflection of the two orders, and the off-diagonal 
elements their coupling. It would seem that in the most general case this matrix is characterized by four 
complex numbers, but by imposing reciprocity this number is limited to four real numbers, and the general 
form of the matrix is 

n^/ ae'^f" ±ji^e'(^o+v-i)/2 \ 

where a, b, and fo^-i are real. Energy conservation provides a relationship between a and 6, but we 
do not need this here. Following the same approach as in the scalar case, we derive in Appendix A the 
transcendental equation for the odd and even fundamental modes. The approach is similar to that of the 
calculation of Lamb waves in acoustics |12j . which also involve two plane waves that couple at the mirrors. 
The resulting equation for triangular lattices reads 

a =F 1 . . 

tani/'oCoti/'-i = (5) 

where ipm = {nkdcos'drn + ¥'m)/2, and the upper and lower signs stand for the even and odd modes 
respectively. For square and rectangular lattices, coti/i-i is replaced by tan^/;_i. 

The next step, which is tedious but straightforward, is to implicitly differentiate Eq. ([s]) to find 

^d^ ^ ^ (ndcos^o + g^) - (ndsec^_i - ndtan^_i sin^p + ^) - ^ff 

'^'^ inrk H^i^^o - - iiirfc Ptan^_i cos^o + l%r)+^dM-o 



3 



where iSq and are related to each other by the grating equation. Note that parameter b does not enter 
this equation. Finally, Eq. ^ may be used in Eq. ^ to express the group velocity in terms of the elements 
of the reflection matrix Q. Even though there are two plane- wave orders, by the grating equation the 
associated f3 = nksiwd differ by a reciprocal lattice vector, and so they lead to the same group index when 
evaluating k^. 

The first terms of Eq. ([6|) in the numerator and the denominator are the same as in Eq. ([3|), the scalar 
case. The second terms are similar, but for the —1 order; as the incident angle of the zeroth order plane 
waves varies, the — f order varies as well, since the angles are related by the grating equation. The third 
terms are new and lead to a phase distortion, rather than an amplitude distortion, of the orders, but we have 
found these to be small in practice. As in the previous section, where there is only one specular reflection, 
Vg = when the denominator in (j6| vanishes. 

Our method can be used to show that Vg — at the BZ edge; we now briefly discuss how. By spatial 
symmetry, the plane wave orders constituting the backwards waveguide mode (denoted by a prime) are 
related to those comprising the forward mode: ??q = — i^q and = — where for consistency we 
have also reversed the grating's period. A' = —A. At the BZ edge, the two grating orders are in Littrow 
configuration, i.e., (3o = tt/A, = — tt/A, and -do = —-d-i. Therefore at the BZ edge, ^'q = 
and — "^o- the backward and forward waveguide modes are composed of the same two plane waves. 
Applying symmetry properties and energy conservation to Eq. Q, we see that the off-diagonal terms of 
R are equal, and that ifo — for rectangular lattices and <y9o ~ + foi' triangular lattices. Also, 
fcot = fc-iti so by solving Eq. (A.3) we find that the amplitudes of the plane wave orders in a mode are 
related: A_i = ±Ao for a square lattice, and — iiAo for a triangular lattice. We have therefore shown 
that two waveguide modes consist of the same two plane waves with the same ratio of amplitudes: thus 
they are indistinguishable. Since by definition the modes propagate in opposite directions, they must have 
w„ =0. 



2.3. Calculation of (po and ip^i 

The waveguide mirrors that we consider are either semi-infinite PCs, or stacks of PC layers in front of 
semi-infinite PCs. The reflection matrix Q describes the amplitudes and phases of the two plane waves' 
reflections and diffractions off this mirror. The matrix, and thus ipo and if-i, may be calculated with 
electromagnetic scattering software by simulating the mirror with one incident plane wave at a time, and 
measuring the complex amplitudes of the reflected plane waves. We calculate the reflection matrix using 
2x2 PC impedance matrices [13 found by an accurate and efficient multipole method [14]. 



3. Results 

We now apply the method from Sec. |2.2| to design a number of PCWs. The most important, and 
computationally slowest, part of our computation is determining the reflection characteristics of the mirrors. 
Since at the beginning of the calculation the frequency and incident angle are unknown, the parameters a, 
ifo, and if-i from Eq. Q need to be calculated for all incident angles and all frequencies inside the bandgap. 
This also provides a key advantage: once these are known the dispersion relation for a waveguide of any 
width can be calculated with very little additional effort. To illustrate this consider Fig. [2ja), for which 
r = 0.3A, where r is the radius and A is the period, n — 2.86, and light is polarized with the H field out 
of the PC plane. The curve labeled "0.5" in this figure, for example, shows the dispersion relation for the 
even mode of a WO. 5 waveguide, i.e., a waveg uide of width 0.5A\/3/2. The other curves give the dispersion 
relations for waveguide of other widths. Using our method it is very easy to calculate the dispersion relations 
of waveguides of many different widths. This fills the entire bandgap, and the color coding gives the group 
velocity for the unique dispersion relation at each point. White regions correspond to bands, where the PC 
does not act as a mirror for that incident angle. The black line is the light line, whereas the blue line is the 



first Wood anomaly. Below the Wood anomaly the scalar treatment from Sec. 2.1 is sufhcient, whereas above 



the Wood anomaly the vector treatment (Sec. 2.2 1 is required. The figure shows the expected behavior: the 



group velocity approaches zero as (3 approaches the edge of the BZ. 
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Figure 2: Dispersion relations for different PCW designs for varying waveguide thickness. Tfie red curves give the waveguide 
width in units of A\/3/2. (a) H out of plane, n = 2.86 and r = 0.3A; black line: light line; blue line: Wood anomaly, (b) 
similar to (a) but with r = 0.404A in the second row of holes: the dispersion relation has a stationary inflection point for a 
Wl waveguide (width d = A\/3/2). (c) Similar to (a), but r = 0.26A and r = 0.34A in the first and second rows of holes, 
respectively: the dispersion relation has a quartic degenerate band edge for d = 0.76A\/3/2. 

In 2008, White et al. 0^ published a waveguide design with a purely cubic dispersion relation, for a 
Wl waveguide. Near this stationary inflection point the waveguide exhibits two desirable properties: the 
modal group velocity is low and efficient coupling into this slow mode is possible [151 116| . The results for 
this design are shown in Fig. [2jb), which has the same parameters as in Fig. [2]ja), except that for the second 
row of holes r = 0.404A. Our results agree well with those of White et al.: we obtain a similar dispersion 
curve for the Wl waveguide, including the presence of a stationary inflection point. Our figure also shows 
how the dispersion relation evolves with waveguide width. 

The final example is a variation on the previous ones in that the radii of the first and second rows of holes 
were changed to r = 0.26A and r — 0.34A respectively. The resulting dispersion relations are summarized 
m Fig.lJc). The fi gure shows that for narrow waveguides {d < 0.7A-\/3/2) the dispersion relation curves 
down when moving away from the BZ edge, while for wider waveguides d > 0.8A-\/3/2, it initially curves 
up. This implies that for waveguides of intermediate width the curvature close to the BZ edge vanishes. 
In other words, the dispersion relation is locally quartic. This is confirmed by the red dashed curve, which 
gives the dispersion relation for a waveguide of width d — 0.76A\/3/2. Like cubic dispersion relations, the 
evanescent modes associated with these features allow for efhcient coupling into slow- light waveguides [17] . 
This example shows the advantages of being able to scan a part of parameter space with little computational 
effort. 

The treatment here is a two-dimensional one in which the structure is taken to be uniform in the direction 
orthogonal to the periodicity. As mentioned, the design of PCWs often starts with a two-dimensional 
treatment to identify the most promising part of parameter space, which is then followed by a full three- 
dimensional treatment. The black curve in Fig. [2jc) shows the results of a fuUy-vectorial 3D eigenmodes 
solution of Maxwell's equations in a plane wave basis |TS], with a slab thickness of 0.58A and the refractive 
index of silicon (n = 3.5). The presence of the quartic dispersion relation confirms that the two-dimensional 
calculations are an excellent predictor for full results. 

4. Discussion and Conclusions 

We have presented a semi-analytic method for calculating the dispersion relation of PCW modes, based 
on the analogy with conventional planar waveguides. This method is efhcient in that calculating the disper- 
sion for different waveguide widths requires negligible additional calculational effort. The most important 
limitation is that our treatment is a two-dimensional one, so that the results for three-dimensional geome- 
tries follow from an effective-index |10| argument. The reason for the two-dimensional treatment is that, 
even for conventional, non-circular waveguides (i.e., not based on PCs) the calculation of the modes for 
three-dimensional structures is usually fully numerical and is not cast in term of mirror properties. In fact 
the key exception to this is when the effective index method is used! The results in Fig. [2jc) confirm though 
that a two-dimensional treatment gives an excellent indication of the full, three-dimensional properties. 
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Other, more minor approximations are that we have assumed that the refractive index of the waveguide 
material does not depend on frequency. The medium dispersion is easily included but this small effect leads 
only to very small changes in the final results. We have also assumed the presence of two propagating 
diffracted orders: the specular order and one other. Though for the parameters of interest no additional 
propagating orders exist, a more complete theory would need to include evanescent orders. However, at the 
frequencies of operation the evanescent orders tend to decay significantly over the width of the waveguide, 
and, additionally, their amplitudes are in fact very small as well. Indeed, we have found that the vector 
theory from Sec. |2.2| is sufficient for our purposes. This is illustrated by the fact that our method reproduces 
the presence of the dispersion relation which locally is cubic (Fig. [2]Jb)), a delicate feature which is sensitive 
to the parameters and to the use of approximations. 
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Appendix A. Derivation of Eq. ([s]) 



In this Appendix we derive Eq. ( 5 1 for a mode consisting of two coupled plane wave orders, using a 
resonance argument and expression (4) for the reflection off the interface with the photonic crystal. We 
do so following an argument set forth earlier for acoustic waves \TT (see also [5]). Since the geometry is 
symmetric the solutions are either even or odd. Let us take the complex amplitudes of the two orders at the 
left of the waveguide, propagating to the right, to be (A_i, Aq) for the —1 and orders, respectively. They 
then propagate to the right-hand mirror and undergo reflection as described by Eq. Q , before returning to 
the left-hand mirror and reflecting again. For a mode to exist, the phase accumulated in this process needs 
to be such that the amplitudes are (A_i, Aq). We can rewrite this condition in vector form as 



RP- ■RP+ 
where R is given by Eq. Q and for square lattices 

P+ = P- = 

where /co.-it = fccos(i9o,-i) 



A-i 











(Al) 



(A.2) 



lattices, we factorize Eq. (A.l), writing P 



is the transverse component of the wavevector. 



RP 



A-i 



± 



A-i 
Ao 



Since P+ = P for square 



(A.3) 



consistent with the result for acoustic waveguides [T^ . 

A subtlety applies to triangular lattices, which necessitates a modiflcation to P. The reflection matrix R 
is defined with respect to a particular point — if this point is moved along the PC interface, then the phases 
of the off-diagonal elements change because Po ^ /3-i. For triangular lattices, which have non-orthogonal 
lattice vectors, the reflection matrices for the left and right mirrors are equal when they are deflned for 
points vertically separated by a half-integer number of PC periods, e.g. A/2 (Fig. A.3). This vertical shift 
must be included in for Eq. ( A.l ) to apply to triangular lattices: 



gi(fe_it<i±^_iA/2) 







„i(feotd±/3oA/2) 



(A.4) 
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Figure A. 3: The PC mirror's reflection matrix is calculated with respect to the point at the centre of the PC's unit cell's edge. 
This means that a phase shift associated with this vertical translation of it A/2 must be included in Eq. ||A.3|. 



By the grating equation, P = P+e *'^o^. Inserting this result into Eq. (A.l) and factorizing, we obtain 



Eq. (A. 3) for triangular lattice PC mirrors but with 



P = 









oikotd 



(A.5) 



Some straightforward manipulation and an application of the grating equation then leads immediately to 
Eq. ([5]), in which b does not appear. In the presence of only the specular order the vectors and matrices 
become scalars and a = 1, leading immediately to Eq. ([T]). 
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